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Abstract 

We have investigated the effect of thermal fluctuation on the interlayer cou- 
pling of three-dimensional Josephson-junction arrays with anisotropic inter- 
actions as a model of layered high-T c oxides, focusing on non-Ohmic current- 
voltage characteristics. Langevin dynamic simulations were performed for 
various Josephson-coupling anisotropics and for various temperatures in both 
zero and finite magnetic fields. We find that, in the highly anisotropic region, 
the interlayer coupling, which causes the non- vanishing critical current (J c ), 
is suppressed by thermal fluctuations, and tends to push up the Kosterlitz- 
Thouless transition temperature (Tkt) of the two-dimensional (2D) regime. 
In this highly anisotropic region, J c begins to decrease drastically near Tkt- 
This suggests that interlayer coupling is screened by thermally excited 2D 
small vortex-antivortex pairs as well as the intralayer logarithmic part of 
vort ex- ant i vortex interactions. Moreover, weak magnetic fields parallel to the 
c-axis affect only the logarithmic-interaction part between the 2D pancake 
vortices. 
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I. INTRODUCTION 



The coherence length (£ a fc)-scale thermal fluctuation in high-T c oxides has been exten- 
sively discussed in association with the thermal excitation of two-dimensional (2D) vortex- 
antivortex pairs. In the strictly 2D case, it is well known that the logarithmic interaction 
between these vortex-antivortex pairs shows power-law behavior of current-voltage (I-V) 
characteristics, and leads to the Kosterlitz-Thouless (KT) transition.0 Recently, the KT tran- 
sition has been observed even in a one- unit-cell(12A)-thick YBCO thin film.! This clearly 
indicates that superconductivity in high-T c oxides is intrinsically two-dimensional. Under 
weak magnetic fields, power-law behavior and the KT transition have also been observed 
in both I-V characteristics and magnetoresistance.il These observations were also confirmed 
by numerical simulations of 2D Josephson-j unction (J J) arrays.0 

On the other hand, high-T c superconductors in bulk are three-dimensional (3D) weakly 
coupled stacks of superconducting Cu02 planes. The existence of a weak but non-zero inter- 
layer Josephson-coupling between Cu02 layers has been observed.il This coupling could alter 
the logarithmic interaction between 2D pancake vortices,i thus blunting transition, shifting 
the Kosterlitz-Thouless transition temperature Tx^JJi and modifying I-V characteristics.iE 
Moreover, the correlation of 2D vortices in adjacent layers is expected to decrease near 
the transition temperature T cr , and it has been observed that layers are decoupled when 
It is not yet known whether layers are completely decoupled just below T cr . 

In this paper, we systematically investigate the screening effect of the interlayer coupling 
on the electric transport properties in both zero and weak magnetic fields (H || c), and 
discuss the results based on the KT mechanism. 

II. MODEL AND FORMALISM 

We start with a 3D anisotropic Josephson-j unction ( JJ) array regarded as a simple lattice 
version of a type-II superconductor derived from a Lawrence-Doniach type model. 

The 
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Hamiltonian is written as 

H = -J ±c cos(A% - Aij) - J\\ c ]T cos{A6 kl -A kl ), (1) 

<i,j>± c <fc,/>|| c 

where A9ij denotes the superconducting phase difference between nodes i and j. Jj_ c and J\\ c 
are intralayer (inplane) and interlayer Josephson- coupling energies related to the magnetic 
penetration depth A|| Cj ^ c , respectively, such that J\\ Cy ± c = $o 2 a./167r 3 A 2 || Ci ^ c , and also related 
to the bare critical current of each Josephson-j unction through /o||c,_i_c=(2e/ft) J\\ c ,± c - Here a 
and $o are a lattice spacing and flux quantum, respectively. We assume that a is comparable 
to the inplane coherence length £ ab . Bare anisotropy ratio 70 = j(T = 0) at T=0 is defined 
as 

70 = 7 (T = 0) = J,,,/ J ±c = (A ±C /A|| C ) 2 (2) 

At finite temperature, 70 can be replaced by effective anisotropy ratio 7 (T). Magnetic 
fields are introduced as frustration 2irf = + Ajk + Am + An around each plaquette with 
A^ = (2e/ftc)// A-dl and / = ifa 2 /$o- The external DC current I is injected by adding 
— (h/2e)IJ2<i,j> ll jA9 i j to Eq. (|l|). We directly calculated the I-V characteristics of the 8 x 
8 x 3 JJ array for various reduced temperatures kBT/J± c , and for various weak magnetic 
fields, via Langevin simulation technique, as previously reported in RefjlO]. 

III. VORTEX DYNAMICS 

Thermal excitation of 2D vortex-antivortex pairs plays an important role in 2D super- 
conductivity (70 = 0). Such a 2D vortex pair separated by distance r within the same layer 
interacts with a logarithmic potential U2D — E± ln(r/£ a &), and leads to the simple power-law 
I-V characteristics: 

R = V/I oc I a -\ (3) 

Moreover, in this 2D regime, the KT transition may be caused by spontaneous unbinding of 
vortex pairs due to thermal fluctuation. At the KT transition temperature T^t, the critical 
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singularity can be observed as a discontinuous jump (called universal jump) from 2 to 
in the temperature dependence of the power- law exponent a — 1. This jump separates the 
non-Ohmic resistive state (a — 1 > 2) for T < Tkt from the Ohmic state (a — 1 = 0) for T > 
T^yjil In the presence of a weak interlayer Josephson-coupling (70 7^ 0), it should be noted 
that the interlayer coupling modifies the logarithmic interaction between vortices within the 
same layer.! The modified vortex-antivortex interaction energy U mo d(r) is distinguished by 
the characteristic size of vortex pairs r 7 (T), which is defined by effective anisotropy ratio 7 
based on Eq. @: 



r 7 (T) = U/VtCO- ( 4 ) 

Here the strong anisotropy regime is defined by r 7 ^> 2a, i.e., y/j <C 0.5, where the thermally 
excited 2D vortex fluctuations are dominant rather than the 3D vortex loop fluctuations.0 
In this regime (^/j <C 0.5), while on the one hand the size r of vortex pairs smaller than 
r 7 adds the r-square interaction to the logarithmic interaction between vorticesjl the larger 
size (r 3> r 7 ) adds the r-linear interaction.^ 

U mod (r) =E C + E 1 (T) ln(r/£ o6 ) + (l/A) 1 (T)E 1 (T)(r/U? (r « r 7 ) (5) 



= E C + E 1 (T) \n(r/U) + (tt/\/2) y j(T)Ei(T)(r/£ ab - 1) (r » r 7 ) (6) 

Here E c is the energy to create a 2D vortex in one layer, and Ei(0) = 2nJj_ c . Compared 
with the second logarithmic term and the third r-square term in Eq. ([5]), the logarithmic 
term is considerably larger than the r-square term, and thus 2D behavior is displayed for 
the smaller-r regime (r r 7 ). For the larger-r regime (r ^> r 7 ), on the other hand, the 
r-linear term is of the same order as the logarithmic term in Eq. (^), and we should expect 
the interlayer Josephson coupling effect to appear clearly in the electric properties. 

Considered with the r-linear form of the Lorentz force, this added r-linear interac- 
tion, based on large-size vortex pair excitations, modifies the I-V characteristics with non- 
vanishing critical current I c in the larger-r regime,^ and this is supported by numerical 
simulations of 3D anisotropic JJ arrays.0 



i?oc {I-I c ) a -\ (7) 

Here, 

/ c ~ yJl(T)cEi(T)/$ U (8) 

and 

a- 1 ~ E x {T)/k B T (9) 

are required. In strongly anisotropic BSCCO with 7 ~ 4 x 10~ 4 , r 7 ~ 1000 A is significantly 
large, and almost pure 2D behavior is expected to appear. On the other hand, in less 
anisotropic YBCO with 7 ~ 4 x 10~ 2 , r 7 ~ 100 A, and this much smaller value of r 7 is 
considered to greatly enhance the r-linear term, even for the 10-unit-cell(120 A)-thick film. 

Focusing on this larger-r regime, we below studied the effect of the interlayer coupling 
on the electric properties in zero and finite magnetic fields for various finite temperatures. 
In particular, we will discuss the power-law I-V characteristics with a non-vanishing critical 
current I c and the power-law exponent a — 1 in Eqs. which eventually result in 

^7(T) and Ei(T) renormalized with the screening effect of the interlayer coupling by 2D 
vortex-antivortex fluctuations in the weak coupling regime. 

IV. RESULT AND DISCUSSION 

The inset of Fig. [I] shows log-log plots of I-R (R = V/I) curves for various reduced 
temperatures k,BT/J±_ c at H = and 70 = 0.04. On the low temperature side (/cbT/Jj_ c < 
0.7), I-R curves display a negative curvature, and seem to obey the modified power-law 
with a non- vanishing critical current I c in Eq. (0), whereas, in the higher temperature 
regime {ksT / J± c > 0.7), they tend to accord with the simple power-law relationship in Eq. 
(0). Also, in this high temperature regime the power-law exponent a — 1 seems to be smaller 
than 2. We can also see this feature more clearly by replotting these curves together with 
Eq. (|7|) as depicted in Fig. 0. Here the dotted lines are fitted for each I-R curve using on 
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X 2 fit based on Eq. (|7|). The I-R curves seem classified into two groups. One group, for 
ksT/J^c < 0.7, has non-Ohmic power-law behavior with a non-vanishing critical current 
I c and with an exponent which a — 1 larger than 2. The other is for kBT/J± c > 0.7 and 
displays nearly Ohmic power-law behavior for a — 1 < 2. 

In the presence of external current, the modified interaction energy between two vortices 
from Eq. is written by 

U mod (r) = E C + E 1 hx{r/U) + E 2 (r/U ~ 1) - (%/c)Ir (10) 
~ E c + E 1 ln(r/U) ~ ($o/c)(/ - h)r, (11) 

where E 2 = */7-^i- F° r the strictly 2D case at 7(T) = 0, the third term in Eq. ( |I0D vanishes 
because E 2 = 0, and Eq. ( |i0|) results in the pure 2D interaction of the logarithmic form.0 
Due to the third term in Eq. (|TTD for 7o 7^ 0, current-induced unbinding of 2D vortex 
pairs are suppressed by 

I c = E 2 (c/U%) = VlMc/^oU)- (12) 

The maximum value of U max and the maximum size r max (U max = U mo d(r max )), that give 
the thermal activation rate T from bound pairs to free vortices by T oc e~ Umax / kBT , become 
larger than in the 2D case. Therefore the resistivity is suppressed in the presence of the 
interlayer coupling. 

Figure § (a) shows the temperature dependence of I c for various values of bare anisotropy 
ratio 70 at H=0. I c decreases with increase in temperature. In particular, I c vanishes in the 
high temperature region for the extremely weak coupling (0 < 70 < 0.1). The effect of the 
interlayer coupling is suppressed by thermal fluctuations, and vanishes to zero. When we 
compare the I-R curves and the temperature dependence curve of I c for 70 = 0.04, we can 
see that Ohmic characteristics with a — 1 < 2 appear in the region where I c vanishes. There 
are two different ranges distinguished by temperature Tj c with I C (T > Tj c ) = 0. Moreover 
I c is extremely suppressed just below T Ic . These regions are believed to be associated with 
the KT transition, as we will discuss later. 
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In the temperature dependence of a — 1 as shown in Fig. |2| (b), we can see a kink-like 
anomaly in the curve for the 2D case (70 = 0) around a — 1 = 2. A discussion of this anomaly 
was based on the KT vortex-antivortex unbinding transition corresponding to a sign for the 
universal jump. For the modest-size (8x8x3) lattice, unfortunately, numerical evidence for 
the KT singularity is inconclusive, since the universal jump is expected to broaden largely, 
according to AT oc 1/ ln(L) with linear size L. From the numerical data, we extract the KT 
transition temperature Tkt 2D , tentatively determined by the exponent oi{Tkt 2D ) — 1 = 2 
for the 2D case at 70 = 0. Even for the weak coupling case (0 7^ 70 < 0.1) we also see 
the analogous kink-like anomalies around a — 1 = 2. This suggests that the KT transition 
would survive in the presence of such a weak interlayer coupling. The temperature at the 
anomaly Txt 3D (1o 7^ 0) seems to increase monotonically with increasing 70 . Moreover, each 
T KT 3D corresponds to each vanishing temperature Tj c in the temperature dependence of I c 
(Fig. (a)), that is, the KT anomaly appears at the time I c vanishes (Tkt 3D (70) = Tc(7o) 
for 70 7^ 0). Since / C (T) oc E 2 (T), both the intralayer logarithmic interaction part Ei(T) 
and the interlayer r-linear part E 2 (T) of Eq. ( |T0"D must be screened by 2D vortex-antivortex 
fluctuations. This screening effect of vanishing to zero of E 2 (T) above Tkt 3D is consistent 
with Monte Carlo results by Minnhagen et. al.0 It is suggested that the system is anisotropic 
3D below T KT 3D and is 2D just above T KT 3D {E 2 {T) ^ for T < T KT 3D , and E 2 (T) = 
for T > T KT 3D ). 

Note that there are plateau regions in the temperature dependence of a — 1 just below 
T KT 3D for the case with the non-zero interlayer coupling. The size of the plateau becomes 
larger with the increase of the interlayer coupling. The plateau begins at the KT anomaly 
temperature Tkt 2D without the interlayer coupling around a — 1 = 2. Just above that 
temperature Tkt 2D , Ic begins to decrease drastically with increase in temperature. This 
rapid decrease in I c can be regarded as the drastic suppression of effective anisotropy ratio 



7(T). This is confirmed by seeing the temperature dependence of yj(T) given by the ratio 
of the exponent a — 1 to the critical current I c : 
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E 2 (T) _ %U I c 
£i(T) ck B (a- 1)T' 



(13) 



Here we use Eq. ([D and Eq. (|12|). Figure |3| shows the temperature dependence of our 



extracted value of yl(T) in Eq. ( |T3"D at = for various bare anisotropy ratios 7 . We 

can see temperature independent behavior in the low temperature phase (T < T KT 2D ), 
i.e., 7(T) has a constant value below Txt 2D ■ This implies that below Tkt 2D both E 2 {T) 
and Ei(T) are renormalized with the same form, and effective anisotropy ratio 7(T) is not 
screened (7(T) = 70 for T < Tkt 2D )- l{T) begins to decrease toward zero above Tkt 2D ■ The 
rapid decrease of 7(T) in the high temperature region agrees with Pierson's Renormalization 
Group (RG) study§ including the screening due to small intralayer 2D vortex pairs. Small 
2D vortex pairs in each plane are also relevant to the weakening of the r-linear part of the 
interaction in the high temperature phase, where the logarithmic part of the interaction is 
fully renormalized. This also indicates that correlation along the c-axis is reduced. 

The qualitative change around 7Vt 3D is interpreted as a sort of dimensional crossover 
from anisotropic 3D vortex fluctuations caused by the Josephson coupling between layers to 
pure 2D vortex fluctuations of the decoupled regime. The decoupled regime above T KT 3D 
exhibits Ohmic-behavior, while non-Ohmic power-law I-R curves of R cx (J — I c ) a_1 appear 
in the weak coupling regime below Tkt 3D - Due to the modest-size (8x8x3) lattice, our 
calculated power-law exponent a — 1 is not exactly 0, but a — 1 < 2, in the decoupled 
regime. Nevertheless here it is interesting that smaller vortex pairs reduced both interlayer 
r-linear and intralayer logarithmic interactions between 2D vortex pairs; full renormalization 
of the logarithmic part of interaction seems to be a trigger for reduction of the r-linear part 
of the interaction. 

We can also see the power-law behavior of I-R curves, fitted with Eq. (|7|), even in the 
weak magnetic field at k B T / J± c = 0.6 and Ha 2 / $0 = 0.0104, in Fig. |4] (a). Log-plus r-linear 
interaction between 2D vortex pairs survives even in the magnetic field. Moreover, there is 
other power-law behavior in magnetoresistance, as depicted in Fig. [| (b), at k B T / J± c = 0.6 
and I/I ±c = 0.4. 
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R oc H 13 (14) 

The coexistence of two power-law behaviors (Eq. (0) and Eq. ([14])) leads mathematically 
to the logarithmic dependence of the power-law exponent a — 1 on the magnetic fields as 
in Fig. |5|: a — 1 oc lnl/i/.i'i The external field suppresses the exponent a — 1 through the 
logarithmic function of the magnetic field Ha 2 /&Q. The focusing points appeared in both 
magnetoresistance and the magnetic field dependence of the exponent a — 1. These focusing 
points have the almost same value of ifa 2 /$o of the order of 0.015. With our assumption 
that a ~ £ o6 , £ ab = 12 A in YBCO gives H ~ 20 T, and £ o6 = 38 A in BSCCO gives 
H ~ 2 T. This suggests that the power-law behavior in YBCO is kept up under ten times 
higher magnetic fields than that in BSCCO. This is consistent with experimental results.0 
In the weak-field regime ifa 2 /$o < 0.015 (Fig. f| (b) and Fig. [^), the picture in terms of 
field-induced unbinding of 2D vortex pairs is relevant even in the strongly anisotropic system 
(\/7o < 0.5). The inset of Fig. || shows effective anisotropy ratio \Jl{T) at k B T/ J± c =0.6 
as a function of Ha 2 /§o. As the magnetic field increases, square-root effective anisotropy 
ratio J j(T, H) ~ E 2 (T, H)/E\(T, H) seems to increase slightly, at k B T / J± c =0.6. From the 
relation a — 1 = Ei(T, H)/k B T in Eq. (||), it is suggested that the external field suppresses 
only the logarithmic part of interaction Ei(T,H), and maintains or enhances the r-linear 
part due to the interlayer coupling E 2 (T, H) ~ dl(T, H)E\(T, H). We can also see this 
tendency in Fig || where I c oc E 2 certainly obeys the \Jj(T, independence, and Ei(T, H) 
(the slope of the curve) screened by the magnetic field is expected to give a larger value of 
^7(T, H) ~ E 2 (T, H)/Ei(T, H) than that in the absence of the field. As a result, effective 
anisotropy ratio 7(T) seems to be enhanced in the 7-constant phase (T < T^t 30 (70, H)), 
that is, the weak magnetic field would apparently strengthen the interlayer coupling. From 
this larger value of 7(T), we expect to find that the plateau in the temperature dependence 
of a — 1, where I c drastically decreases, is larger under the weak magnetic field than in the 
absence of any field. Because of this clearer plateau, the universal jump could survive in the 
magnetic field, as observed experimentally! 
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In summary, we have performed Langevin simulations of the 3D anisotropic JJ arrays 
for a model of layered high-T c oxides. We found that the effect of vortex fluctuation due 
to small inplane vortex pairs effectively weakens the interlayer coupling and wipes it out, 
that is, effective anisotropy ratio 7(T) is screened by thermally excited 2D vortex-antivortex 
pairs in each layer. Just above Tkt 3D , layers are decoupled, and the system becomes 2D. 
This is consistent with Monte Carlo study by Minnhagen and Olsson0 and with RG study 
by Pierson0. Signs of the transition seem to appear around the temperature Tj c where the 
critical current I c vanishes. Moreover, we can see that Tkt 3D (70) tends to shift slightly up in 
the high temperature direction for small 70 7^ 0. When the interlayer coupling becomes small 
for < 70 < 0.1, I c vanishes in the high temperature region, and this vanishing of I c may 
play a role in triggering the KT transitions. Since 7(T) drastically decreases with increase in 
temperature just below Tkt 3D , such a weakening effect of the interlayer coupling must mainly 
be due to screening by the high density of 2D vortex pairs. Thus, at around Tkt 3D = Tj B , 
there appears a sort of dimensional crossover from the 3D anisotropic system to the pure 2D 
system. Magnetic fields affect only the logarithmic interaction between 2D vortex-antivortex 
pairs. In our simulation, we have not introduced the effect of the coherence length £ c . Near 
T c , significant growth of £ c beyond the spacing between Cu02 layers is expected to give 
less anisotropic or isotropic 3D behavior far above Tkt 3D for \Jl/(l — T/T c ) ^> 1. The 
KT transition observed in YBCO appears to be a sort of dimensional crossover from the 
anisotropic 3D system to the strict 2D system, while the one observed in BSCCO is the 
typical KT vortex-antivortex unbinding transition in the pure 2D system. At that point, 
YBCO in bulk can be distinguished from BSCCO. 
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FIGURES 

FIG. 1. Resistance R = V/I as a function of / — I c at 70 = 0.04 for various values of reduced 
temperature kBT/J± c = 0.6 in the absence of a magnetic field. Inset: resistance R as a function 
of external current /. 

FIG. 2. Temperature-dependence for various values of bare anisotropy ratio 70: (a) Critical 
current I c as a function of temperature kBT/J± c , (b) Power-law exponent a — 1 as a function of 
temperature &bT/Jl c . 

FIG. 3. Temperature dependence of E2 (T) jE\ (T) for various values of bare anisotropy ratio 
To- (E2 (T)/Ei (T) is proportional to the square-root of j(T)) 

FIG. 4. Power-law type behavior for various values of bare anisotropy ratio 70 at 
kBT/J± c = 0.6: (a) Power-law magnetoresistance at current I/Iq = 0.4, (b) Power-law type 
I-R curves at magnetic field Ha 2 /&o = 0.0104. 

FIG. 5. Logarithmic behavior of the power-law exponent a — 1 at ksT/ J± c = 0.6;the exponent 
a — 1 as a function of magnetic field Ha? Inset: E2/E1 as a function of Ha 2 /^Q. 

FIG. 6. Anisotropy dependence at /cbT '/ 'J j_ c = 0.6 for various magnetic fields; critical current 
I c as a function of VtOH with effective anisotropy j(T). 



13 



